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1. INTRODUCTION 

The interest in single-stage topologies in the last years has been grown due to the need to reduce 
size, cost and increase the efficiency of the power converters [1], [2]. The single-stage topology must 
perform two functions: 1) controlling the dc voltage and 2) ensuring that the output AC voltage is sinusoidal 
with low THD. Z-source inverter (ZSI) is a single-stage converter with buck-boost characteristic property. 
All switching states of the VSI is used in the ZSI in addition to a shoot-through state. ZSI suffers from many 
drawbacks such as, unidirectional power flow, high inrush current during starting, and discontinuous input 
current [3]. The researchers try to solve these drawbacks in two ways. The first way, by changing the 
configuration of the impedance network [3]-[6]. The second was, by improving a proper control strategy for 
various application requirements [7]. The number of passive components in the impedance network of ZSI 
such as inductors and capacitors are relatively large. 

Split source inverter (SSI) is a single-stage boosting converter which has different circuit topologies 
in literature [8]-[12]. SSI has the capability of using the same switching states as the VSI for its basic 
operation and there is no need for extra states [8]. SSI provide advantages over the existing single-stage 
topologies as: 1) continuous dc-link voltage, 2) continuous input dc current, 3) lower passive components 
count, 4) using the same standard modulation schemes as the VSI for its basic operation. On the other hand, it 
suffers from the following disadvantages: 1) a sufficient dead band time should be generated, 2) unequal 
current distribution among the different switches [9]. 
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In literature, SSI uses traditional PWM with PI controller which suffers from the low response and 
minimum phase [9]-[12]. Model predictive control (MPC) is presented with the high-performance response 
and no minimum phase problem [13]-[15]. MPC is a controller, which uses the discrete-time model of the 
system to predict the future behavior of the controlled variables for all possible voltage vectors of the 
converter. The predicted variables are then used to obtain optimal action by choosing the voltage vector that 
minimizes a predefined cost function [16]. Hence, it is more suitable for the discrete nature of the power 
electronics converters, and it lends itself easier to its application [17]-[20]. MPC has many advantages such 
as: 1) it is easy to understand, 2) it can be applied to a wide variety of systems, 3) it deals with multivariable 
cases, and nonlinearities [21], and 4) capable of using constraints [22]-[24]. 

In this paper, finite control-set model predictive control (FCS-MPC) is modified to control the AC 
and DC sides of the single-phase SSI. The cost function consists of output AC voltage error, DC-link voltage 
error, and DC input current error. Hence, only one cost function has been used to control the three variables. 
For proper operation of the modified MPC, weighting factors are used within the cost function. Harris Hawks 
optimization (HHO) technique is used to determine the optimal values of the weighting factors. 


2. FINITE CONTROL-SET MODEL PREDICTIVE CONTROL FOR SPLIT-SOURCE INVERTER 

The single-phase SSI employed in this paper is shown in Figure 1 [10]. The principle of operation is 
described with the aid of Figure 2 and Figure 3 as follows. When the switch Qs is turned ON the inductor Lac 
is charged (see Figure 2), and the SSI can operate with all possible output voltage (+Viny, -Vinv, and zero). The 
inductor Lac is discharging into the capacitor (Cac) when the switch Qs is turned OFF. At this moment, the 
SSI can operate with only two possible outputs (+Vin, and zero) as shown in Figure 3. The DC-link voltage 
(Vin) is function of the duty cycle of switch Qs and is given by (1). 


Va 
Vinv = ma (1) 


In order to design FCS-MPC, the following steps are required: 1) the system model, 2) the discrete time 
model, and 3) the cost function. 


Figure 2. Possible switching states of the single-phase SSI at charging Lac 
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Figure 3. Possible switching states of the single-phase SSI at discharging Lac 


2.1. System model and discrete time model 
2.1.1. AC side model 
The governing differential equations of the LC filter variables are: 


di 

Lyf = v4 -% (2) 
dve . . 

Cp TE = ip — tg (3) 


where, ve is the AC load voltage, v; is the inverter output voltage which applied to the filter, i is the filter 
inductor current, i, is the AC load current, Cy is the filter capacitance, and Lyis the filter inductance. 
The equations can be rewritten in a state-space system form as: 


dx 


q Ax + Bu; + Baio (4) 
where, 
x= [ly UV]? (5) 
0 = 1 0 
Lf a E 
A=|, B= |. and By = |= (6) 
— 0 0 Cr 
Cf 


The output voltage expressed in a vector-matrix form is given as: 

ve =[0 1]x (7) 
The discrete time model is obtained from (4) for a sampling time Ts and expressed as: 

x(k + 1) = Agx(k) + Byvi(k) + Bagio(k) (8) 


where, 
Ay = e*%, By = J, e*Bdt, and Baq = f° e“ Badt (9) 


2.1.2. DC side model 

The DC side of the single-phase SSI consists of an input inductor (Lac) and DC-link capacitor (Cac). 
The governing equations of the DC side has two cases depending on the state of switch Qs [25]-[27]. Case 
one, the switch Qs is turned on and the input inductor (Lac) is charging as shown in Figure 4(a). The dynamic 
equation of the input inductor current and the DC-link capacitor voltage can be presented as: 


di c 
Vac = Lac £ (10) 
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aV inv 


ip = Cac i (11) 
Where, Vac is the input DC voltage, iac is the input DC current, and Vin is the DC-link voltage. 
For first order system, the discrete time model is obtained by simple Euler approximation as: 
pass see ea (12) 
dt Ts 
Hence, the discrete time model of the input inductor current and the DC-link capacitor voltage are: 
; ; Ts 
bac(k + 1) = iac (k) + 7 Vac (k) (13) 
Tsis 
Vinv (K + 1) = Vim (k) + Tac f (14) 


Case two, the switch Qs is turned off and the input inductor (Lac) is discharging as shown in Figure 4(b). The 
dynamic equation of the input inductor current and DC-link capacitor voltage can be presented as: 


Vac — Vinv = Lac £ (15) 

iac — ip = Cac = (16) 
The discrete time model of the input inductor current and the DC-link capacitor voltage are: 

iac(k + 1) = iac (k) + 7 > Vac (k) — Vinw KD) (17) 

Vinw(K + 1) = Vins (k) + <= (iac (©) — ip(k)) (18) 


Figure 4. The two cases of the SSI circuits at (a) Qs is turned on and (b) Qs is turned off 


2.2. Model predictive control strategy 

In this section, the MPC controller is modified to generate the switching action for the single-phase 
SSI switches according to input power, DC-link voltage, and AC output voltage. The duty cycle of switch Qs 
is determined inside the MPC controller. The modified MPC algorithm for the single-phase SSI system is 
shown in Figure 5. At any instance (k) and from knowing the current states of the DC-link voltage (Vin) and 
DC input current (iac), these variables can be predicted in both ON state and OFF state of the switch Qs. In 
the ON state, The DC input current (iac) and the DC-link voltage (Vin) are predicted from (13) and (14). The 
single-phase SSI can operate with all possible voltage vectors as shown in Figure 2. The cost function gon is 
used to calculate the error between each of the predicted variables and their reference values. The cost 
function consists of the AC side error, the DC-link error, and the DC input current error as in (19). Weighting 
factors are used for the proper operation of the modified MPC. Harris Hawks optimization (HHO) technique 
is used to calculate the weighting factors. The voltage vector that minimized the cost function (gon) is stored. 
In the OFF state, the DC input current (iac) and the DC-link voltage (Vin,) are predicted from (17) and (18). 
The single-phase SSI can operate with two possible voltage vectors as shown in Figure 3. The cost function 
(Zo) is used to obtain the optimal voltage vector at OFF state as in (20). In both ON and OFF states the 
output AC voltage is predicted from (8). The optimal voltage vector from the ON state is compared with the 
optimal vector of the OFF state, and the vector that has minimal error is chosen to apply in the next sample. 
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Jon = w, |v", T; ve(k + 1)| t W2|Vinv i Viny (kK + 1)| as Wsliac” = lac (k + 1)| (19) 
Ioff = w,|v", a Ue (k T 1)| T wz |Vinv x inv (k + 1)| + wsliac = lac(k F 1)| (20) 


where, w; is the weighting factor of the AC load voltage error, w2 is the weighting factor of the DC-link 
voltage error, and w3 is the weighting factor of the DC input current error. 


k1 Measure x(k), Vav- Ide 
Estimate ref .signals 
Predict tg Vnv Predict Lac, Vinv 
From (17), (18) From (13), (14) 
Lon(®, Vinstac) Zon (X, Vavta 
From (20) From (19) 
Optimization Optimization 
S opt OFF Sopt_ON 
Yes ` 
No 
Sopt(k+ 1)=Sopt_orr Sopt(k+ 1)=Sopt_on 
k+l Store Sopt(k+1) 


Figure 5. FCS-MPC algorithm 


3. HARRIS HAWK OPTIMIZATION 

Harris Hawks optimization (HHO) technique is used to obtain the weighting factors of the cost 
function. The integration of absolute error (IAE) of AC load voltage error, evac, the DC-link voltage error, e€vac 
and the error current in the DC input current, eiac, are used as the objective function (21). As shown in Figure 6, 
HHO mimics the Harris hawks’ cooperation behavior in chasing the prey, which is called surprise pounce [28]. 
The optimized values of the weighting factor are the prey that hawks tend to catch. As any optimization 
technique, HHO consists of the exploration and exploitation phases. In the exploration phase, the hawks intend 
to detect and track the prey by searching in some random areas. As the prey energy decreases, the algorithm 
transits from exploration phase to exploitation phase. In the exploitation phase, the Hawks attack the intended 
prey detected in the previous phase by performing the surprise pounce [29]. 


IAE = Jy levac(t)|at + Jy \evac(t)|at + Jy leiac(t) at (21) 
where, 
epac = Ve” — Ves Cvac = Vinu. —Viny, and eiac = base — lade (22) 
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INPUT Number of search agents 7, the number of iterations 7, dimensional | Where 
space D Es is the initial state of the prey energy 
Initialize the position of hawks F; (i=1, 2..., 2) . : £ 
WHILE the end condition is not satisfied is the current iteration counter 
FOR i=1:n is the max iteration 


aaa the objective function value ofeach hawk through is a randoni pmaber iinde (1); 


ins P = a indicates the position of the hawk 
Sa gaat p the position of the intended prey Forays according to the position of other 


FOR iaia team members 
Update the initial energyZo, jump strength J Cs, Cz, C3, Ca are random numbers inside 
Update the escaping energy (0,1) updated in each iteration 
E=2E a-b) YW is the current position vector of 
T hawks 
n - s 
IF E |> 1 (Exp loration phase) Yana(t) is a randomly selected hawk from the 
IF K=randQ<0.5 current population 
Update the position of the hawk ¥.(0 is the av -i f the hawks 
Y ¢+D=0 . (3 .@))-c,(LB +c,(UB -1B)) =(t) is the average position of the hawks, 


ELSE LB &UBae the lower and upper bounds of 
Update the position of the hawk variables 
¥ ¢+D=!__@-¢. F , @)-2eF @| is the chance of a prey in successfully 
END IF escaping 
END IF 


IF jE \<1 (Exploitation phase) . ; 
IF 20.5 and E |205 (Soft besiege) is a random vector by size 1xD 
720.5 an 20.5 oft besiege, : ` ; ; 
the Lévy flight functi 
Update the position of the hawk ee a 
T @+D=@ 0) -E PE- ©) tanks a E escape. 
ELSEIF 205 and |E|<0.5 (Hard besiege) Ba S 
Update the position of the hawk 
F ¢+D=!..@-EF.0O7T O 
ELSEIF r<0.5 and |E |20.5 (Soft besiege 


progressive rapid dives) 
Update the position of the hawk 
beni f E) -E PE me GT @)| f F(H)<F(T (t)) 
| H+SxLF(D) if F(G)<F(F (t)) 
ELSEIF 7 <0.5 and E|<0.5 (Hard besiege with 
progressive rapid dives) 
Update the position of the hawk 
r PETS w) -E PE mlt) -E nO f F(H)<FE (t)) 
j H+S*xLF(D) if F(G)<F{¥ (t)) 


is the dimension of the problem 


END IF 
END IF 
Check and correct the new positions based on the boundaries of variables 
END FOR 
END WHILE 
RETURN Fores, Which represents the optimal values 


Figure 6. Pseudo code of Harris Hawks optimization 


4. SIMULATION RESULTS 

The single-phase SSI, shown in Figure 1, is simulated using MATLAB/Simulink. The system 
parameters shown in Table 1. The modified MPC is used to control the single-phase SSI where all the five 
switches (Q1-Q5) are controlled. The AC output reference voltage is chosen to be 220 V, 50 Hz. The DC- 
link reference voltage and the input reference current are chosen to be 500 V, and 8.5 A (1200 W) 
respectively. The weighting factors (in (19) and (20)) are obtained by HHO and equal to w;=1, w2=0.16, and 
w3=0.11. The AC output voltage waveform (vc) is shown in Figure 7. The output voltages are sinusoidal with 
low THD of 1.49% and peak value equal to 311 V as shown in Figure 8. Figure 9 shows the load currents 
(io), which equals 5.5 A rms. The load consumed nearly 1200 W from the supply. The DC-link voltage (Vim) 
is shown in Figure 10, and equal to 495 V with low ripples (less than 2%). The DC input inductor current 
(iac) is Shown in Figure 11 and has average value equals 8.5 A. The results show that more than one variable 
can be controlled with MPC with proper selection of weighting factors. And this is the most important feature 
of the MPC over other controllers. 

The dynamic behavior of the single-phase system is investigated by applying sudden change in the 
load from 40 Q to 20 Q (100% load change) at 11.05 sec as shown in Figure 12. The figure shows that the 
current undergoes increase from 5.5 A to 11 A (rms) and take 0.1 msec to reach its steady-state value. The 
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load power increases from 1200 W to 2380 W. The output voltage (ve) waveform is shown in Figure 13. The 
figure shows that there is a slight effect on the load voltage due to this sudden load change. The THD of the 
output voltages is 1.55% as shown in Figure 14. The DC-link capacitor voltage (vim) is shown in Figure 15. 
The figure shows that the DC-link voltage is slight affected by load change. The boosting inductor current 
(iac) is Shown in Figure 16, which its average value has increased from 8.5 A to 17 A. 


Table 1. Single-phase SSI system parameters 


AC side 

Reference output voltage ve” 220 V rms 

Frequency f 50 Hz 

Switching Frequency fe 10 KHz 

Output filter inductor Ly 5 mH 

Output filter capacitor Cp 100 uF 

Load resistance R; 40Q 

Load inductance L; 10 mH 

DC-Link side 

DC-source 140 V 

DC-link voltage 500 V 

DC-link capacitor Cy. 2000 uF 

DC-link input inductor La. 3 mH 
400 - 3 
200 J7 

č o | 
om 
-200 
-400 - L 1 1 I d 
10 10.05 10.1 10.15 10.2 10.25 10.3 10.35 10.4 
time (sec) 


Figure 7. The AC load voltage (ve) 


Selected signal: 600 cycles. FFT window (in red): 20 cycles 
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Figure 8. THD of the AC load voltage (ve) 
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Figure 9. The AC load current (io) 
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Figure 10. The DC-link voltage (Vim) 
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Figure 11. The DC input current (iac) 
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Figure 12. The AC load current (io) at 100% step load change 
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Figure 13. The AC load voltage (vc) at 100% step load change 
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Figure 14. THD of the AC load voltage (vc) at 100% step load change 
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Figure 15. The DC-link voltage (Vim) at 100% step load change 
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Figure 16. The DC input current (iac) at 100% step load change 


5. CONCLUSION 

In this paper single-phase split-source inverter (SSI) has been controlled by a modified FCS-MPC 
based on Harris Hawks optimization. The AC load voltage, the DC-link voltage, and the DC input current are 
controlled simultaneously. Weighting factors determine the importance of each term within the cost function 
and are obtained using Harris Hawks optimization technique. The SSI performance is validated under a 
certain load, which consume 1200 W. The results show that the AC load voltage is well regulated with low 
THD (1.49%), and the input power is controlled while maintaining the DC-link voltage content with low 
ripples (2%). The system is also simulated under 100% load change, which the load consumes 2380 W. The 
results show that FCS-MPC has a high-performance response. In both cases, the DC input current has high 
ripples (8.5 A ripples) as it has a low weighting factor. As the DC input current weighting factor increases, as 
the THD of the AC voltage and ripples of the DC-link voltage increases too. The FCS-MPC algorithm 
sacrifices the DC input current to regulate the AC and DC voltages. 
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